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Starting from a continuous time random walk (CTRW) model of particles that may evanesce as 
they walk, our goal is to arrive at macroscopic integro-differential equations for the probability den- 
sity for a particle to be found at point r at time t given that it started its walk from ro at time t = 0. 
The passage from the CTRW to an integro-differential equation is well understood when the parti- 
cles are not evanescent. Depending on the distribution of stepping times and distances, one arrives 
at standard macroscopic equations that may be "normal" (diffusion) or "anomalous" (subdiffusion 
and/or superdiffusion) . The macroscopic description becomes considerably more complicated and 
not particularly intuitive if the particles can die during their walk. While such equations have been 
derived for specific cases, e.g., for location-independent exponential evanescence, we present a more 
general derivation valid under less stringent constraints than those found in the current literature. 

PACS numbers: 02.50.Ey,82.40.-g,82.33.-z,05.90.+m 



X 



I. INTRODUCTION 

Continuous time random walks (CTRW) offer a sweep- 
ing framework to describe the dynamics of particles 
whose motion may be "anomalous." That is, in addi- 
tion to providing a way to describe the motion of diffu- 
sive particles at a more microscopic level than, say, the 
more macroscopic diffusion equation, CTRWs can also 
be used to describe particles whose motion is subdiffusivc 
or superdiffusive. The connection between CTRWs and 
the associated "more macroscopic" description, which is 
in general an integro-differential equation rather than a 
(partial) differential equation as in the case of diffusion, 
has been firmly established. The more macroscopic de- 
scription is often attractive because all the machinery 
of integro-differential equations and associated bound- 
ary conditions can be brought to bear in their solution. 
It should be noted that the starting point for problems 
that involve subdiffusive or superdiffusive particles is of- 
ten simply the integro-differential equation itself. How- 
ever, one must exercise caution in simply starting with 
such a macroscopic description. It should also be noted 
that these integro-differential equations can arise from 
microscopic models other than a CTRW, the latter not 
always being the appropriate framework. An example is 
the motion of particles in random landscapes of various 
sorts, for example, one in which each site is associated 
with a potential well whose depth is chosen from a ran- 
dom distribution. The walk in a landscape of potential 
wells of random depths may be normal in the sense that 
the escape from each well follows a Kramers law, but the 
irregularity of the landscape may give rise to trapping 
events that may slow down the progress of the particles 
to the point of subdiffusion. 

The situation becomes much more complicated when 
the moving particles also undergo reactions [l|, ■ This 
applies to situations such as reversible or irreversible con- 



version to a different species (e.g., A — > B [3, d, Q or 
A ^ B |), reactions giving rise to propagating fronts 
(say, A + B -> C Q, A + B -> 2A Q), binary reac- 
tions (e.g., A + A — > 0), or even spontaneous evanescence 
{A — > 0). Focusing on the latter case, the description of 
the evolution of such evanescent particles at the macro- 
scopic level of a subdiffusion or superdiffusion equation, 
and the investigation of the proper way to include the 
evanescence or reaction in such equations, is a matter of 
continuing study and is usually carried out in the context 
of rather specific models [1, 

In this paper we pursue this goal somewhat more 
generically: we consider the deduction of subdiffusion or 
superdiffusion equations when the moving particles are 
intrinsically evanescent. The death rate of the particles 
may in general depend on location, time, how long a par- 
ticle has spent at a particular location - one can think 
of a large variety of death scenarios. Our presentation 
is based on a CTRW formulation generalized to include 
particle evanescence. One step of the discussion concerns 
the way in which evanescence can be built into a CTRW 
model, a choice which is not unique. The second step 
is then to go from a CTRW model to a subdiffusion or 
superdiffusion equation as appropriate. 

In Sec. [IT] we start with some basic CTRW quantities 
and relations, which we then use in Sec. IHII to construct 
an integral equation for the probability density of finding 
a particle at a certain position at time t given that it be- 
gan its motion at a different position at time t = 0. This 
integral equation includes the possibility that the parti- 
cle evanesces along the way, and is the starting point for 
the further derivation of subdiffusive and superdiffusive 
equations. In Sec. IIVI we explicitly construct our frac- 
tional equations for the case of subdiffusion, superdif- 
fusion, and a mixture of the two. We conclude with a 
summary and some thoughts for the future in Sec. [V] 
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II. STARTING WITH A CTRW 

Our goal is to arrive at integro-differential equations 
for the probability density w(r, t\r , 0) for a particle that 
started its journey at point r at time t = to be at r 
at time t. No matter how we arrive at such an equation, 
the description of the motion of the particles as a CTRW 
requires the introduction of the probability density ^(r — 
r', t—t') that a random walker jumps in a single step from 
r' to r after waiting a time interval t—t' at the position r'. 
The dependence on only the difference of the two position 
vectors reflects an assumption of spatial homogeneity for 
the jumping mechanism. Associated with this probability 
density is the usual waiting time probability density, 

V(*) = J dr'*(v',t), (1) 

and also the probability density for a single step displace- 
ment, 

/>oo 

X (r)= / df*(r,f). (2) 
Jo 

The "normal" or "anomalous" (subdiffusive or superdif- 
fusive) character of the process depends on the forms of 
the probability densities, which we will specify later. 

Since we deal with evanescent particles, we also intro- 
duce the probability E(r, r'; t, t') that the particle has not 
died spontaneously during the stepping process described 
by ^(r — r', t — t'). Also, finding the particle at r at time 
t does not mean that it jumped onto that location ex- 
actly at that time. In fact, it might have jumped there 
at an earlier time and then just waited there without dy- 
ing, or it might have jumped there at an earlier time, 
returned any number of times before i, and then waited 
there. As a starting point in our route to an integral 
equation for w(r, t|ro, 0), it is thus appropriate to intro- 
duce 5n( r ) t|ro, 0), the probability density that a particle 
be at position r at time t exactly after making its nth 
jump, given that it started its walk at time t = at po- 
sition ro 0. This initial condition is described by the 
equation 

(fc(r,i|r 0) 0) = 5(r-ro)«(t). (3) 
From this follows the definition 

oo 

g(r,t|ro,0) = 5^g n (r,t|r ,0), (4) 

n=0 

which is the probability density that the particle steps 
onto position r exactly at time t (regardless of how many 
times the particle has stepped anywhere, including on r, 
before this time). 

It is straightforward to write integral equations for the 
probability densities q and q n . For q n no further dis- 
cussion beyond the definitions is necessary to write the 



integral recurrence equation 

<Z„+i(r,t|r ,0) = / dr' J dt' ¥(r - v't - t') 

x gn (r',i'|r ,0) S(r,r';t,t'). (5) 

In words, this says that the probability that the particle 
arrives at location r at time t on the (n+ I)st step is equal 
to the probability that it arrives at r' at an earlier time 
t' on the nth step and then steps from r' to r at time t, 
provided it does not die in the interval between these two 
steps. The usual integral equation without evanescence 
is immediately recovered upon setting 3 equal to unity. 
Summing this equation over step number n, the integral 
equation 

g(r,t|r ,0) = / dr' J dt' *(r - r', t - t') 

x q{r', t'\r , 0) S(r, r'; t, t') + q (r, t\r 0> 0) 

(6) 

immediately follows. 

The more complex question is now how to go from 
these integral equations to one for the desired probability 
density w(r, t\ro, 0) and thence to a fractional diffusion 
equation. This is our pursuit in the next section. 

III. INTEGRAL EQUATION WITH 
EVANESCENCE 

Before obtaining an integral equation from the CTRW 
setup of the last section, we note that one could choose 
to write such an equation directly, for example, 

w(r,t|r ,0) = J dr' J dt' T(r,r';t,t') 

xw{r',t'\r Q ,0)S(r,r';t,t') 

+ $(t)S(r,r;t,0)5(r-r ), 

(7) 

where we have introduced 

pt poo 

$(t)=l-/ dt'ip(t') = dt>(t'), (8) 
Jo Jt 

the probability that the particle does not take a step in 
the entire time interval up to time t. Note that we have 
made a point of using 5 rather than 5 for the function in- 
dicating that no death occurs, because in general there is 
no reason to expect these two to be the same. The func- 
tion S that appears in the integral equations obtained 
above is the probability that the particle does not die in 
a time interval exactly delineated by two steps (and none 
in between), one taking it to location r' and the other to 
location r. On the other hand, S is the probability that 
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the particle does not die in a time interval t — t' as it 
moves from r' to r, with no reference to steps. There is 
no a priori reason for these two probabilities to be equal. 
If particles die at a constant rate independent of position 
and of when steps take place, then these two probabil- 
ities would be equal. However, for example if particles 
die only when they take a step, or, conversely, if parti- 
cles are more likely to die if they remain at one location, 
then these probabilities would not be equal. Also, we 
have denoted the kernel of the integral equation by the 
new symbol T because a priori we do not know its con- 
nection to ty. One might attempt reasonable guesses, but 
that is all so far and in fact somewhat risky. 

Our goal is to obtain an equation such as Eq. ((7J) start- 
ing from a CTRW (rather than just writing it down). We 
therefore begin with Eq. © together with the exact re- 
lation 

itf(r,t|ro,0)= / dt'$(t-t')q(r,t'\r ,0)S(r,r;t,t'). 
Jo 

(9) 

This exact relation simply says that the probability den- 
sity for the particle to be at r at time t is the probability 
that it stepped onto that location at time t' < t and then 



neither moved nor died until time t. Earlier arrivals and 
returns are implicitly included in this relation. 

In order to proceed toward an integral equation of the 
form ([7]) with known kernel we find ourselves having to 
make two admittedly debatable assumptions. One is that 
the functions S and S are equal. The other is that this 
function can be written as a ratio of the form 



S(r,r';i,*') = 



<p{v>,t>) 



(10) 



The special case of location-independent exponential de- 
cay ip(r, t) = exp(— kt) (that is, the case of evanescence at 
a constant rate, mentioned earlier), is that of Sokolov et 
al. [H, and in this case 3(r, r';t, t') = exp[— k(t — t')]. 
This is the only choice for which S depends on time 
only through the difference of its time arguments. While 
Eq. (flU)) is limiting, it is more general than the cases that 
have been treated in the literature. 

To make use of the special form (ITU|) we start by mak- 
ing the replacements i — ► t' and t' — * t" in Eq. ©. We 
then multiply by $(i — i') E(r, r; t, t'), integrate over t', 
and use relation ([9]) to write 



w(r,t|r o ,0) 



jT dt' $(t - t') S(r, r; t, t') ^ J dr' J* dt"^{r - r', t' - t")q(r' , t"\r , 0) S(r, r'; t', t")j 

f dt' $(i - t') H(r, r; t, t') q (r, t'\r , 0). (11) 
Jo 



I 

This equation can be manipulated through a number of on the right side and divide both sides of the equation 
simple steps. First, we insert the special form Eq. (flQ]) by ip(r,t). This leaves us with the equation 



w(T,t\r ,0) 
<p(r,t) 



, , , n q(r',t" |r ,0) 



Jo 



<p(r,t>) 



(12) 



Next, we recognize that the first term on the right is taking Laplace transforms. Multiplying the rearranged 
a double convolution with respect to time. The mem- equation by tp(v, t) it is then easy to see that 
ory kernels can then be exchanged, as can be verified by 
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tfl(r,t|r o ,0) = J dr'^^'*(r-r',t-0^|y ([' dt" *(t' - t")q{v' , t"\r , 0)^^ 



r 



(13) 



Finally, using Eq. ([5]) and the initial condition (J3J), and sired integral equation 
again recalling the special form (fT0|) we arrive at the de- 



w(r,t|r o ,0) 



dr' / dt' *(r -r',t- t')w{r', t'\r , 0)3(r, r'; i, t') + $(i)£(r, r; t, 0)5(r - r ). 



(14) 



r 



Equation (]14|) is the starting point for the derivation of 
various fractional diffusion equations for different forms 
of the single step probability density ty(r — r', t — t') of 
the underlying CTRW. It is therefore a centerpiece of this 
work. Note that the special form (TIU)) has led to a kernel 
in the integral equation that is precisely this single step 
probability density even though the times t and t' are not 
necessarily associated with jumping times. We also stress 
once again that in addition to the simplification (|10[) we 
have assumed the equality of the survival functions 2 and 
S. Equation (| 14[) has a clear physical interpretation: it 
considers all possible ways for a particle to be at point r 
at time t by looking at the positions r' at prior times t' 
and then tracking their subsequent arrival at the desired 
point. If the point under consideration is the initial posi- 
tion, the equation tracks the possibility that the particle 
has not moved by time t. In the language of Hughes [l3j] 
(Sec. 3.2.8), Eq. (fl"4]) corresponds to a "partition over 
the last step." The equation counts only those particles 
that do not evanesce in the process. 

Finally, we shall implement one additional simplifying 
assumption widely adopted in the literature, namely,that 
the waiting time and jump displacement distributions are 
mutually independent, so we can write 



tf(r,t) = VWxM- 



(15) 



Different fractional diffusion equations then arise depend- 
ing on the behaviors of the tails of these distributions. 
We proceed to present various cases in the next section. 
Our derivations closely follow known results presented in 
a number of helpful review sources such as the reports 
of Metzler and Klafter [1JJ, fijl and a recent multiau- 
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thored compendium on anomalous processes [12 ]. Our 
main purpose here is to add evanescence to the mix and 
to determine how the evanescence "reaction" enters these 
equations. In particular, when tp(t) has long tails and 
x( r ) does not, we will arrive at a fractional subdiffusion 
equation. When, on the other hand, x(r) has long tails 
but ip(t) does not, we arrive at a superdiffusive equation. 
The most "anomalous" case occurs when both have long 
tails, which leads to a bifractional equation. 

We end this section with a practical consideration. In- 
stead of working with the probability density of interest, 
w(r, t\ro, 0), it turns out to be more convenient to work 
with a ratio introduced earlier, 



7?(r,i|r ,0) 



w(r,t\r ,0) 
P(r.«) 



(16) 



Dividing Eq. (I14|) by f(r,t) and taking the Fourier (for 
space)-Laplace (for time) transform, we find 



??(q, u) = x(q)^(w)??(q, u) + $(u) 



Together with the relation 



*(„) = i_m 

u u 



<^(r ,0) 



(17) 



(18) 



and the convolution theorems for both Fourier and 
Laplace transforms, we arrive at the Fourier-Laplace 
transformed reaction-diffusion equation with evanes- 
cence, equivalent to Eq. (| 14[) when the memory kernel 
can be factorized as in Eq. (fT5"]) . 



ufj(q, u) = ux(q)-0(u)?7(q, u) + (l - ^(u)J 



(19) 
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IV. FRACTIONAL DIFFUSION EQUATIONS 
WITH EVANESCENCE 



In this section we proceed to deduce the fractional dif- 
fusion equations appropriate for long-tailed waiting time 
distributions, for long-tailed jump distance distributions, 
and for both simultaneously. 



A. Fractional subdiffusive equation with 
evanescence 

Subdiffusion is characterized by a waiting time distri- 
bution with a tail so long that it lacks integer moments, 
that is, 



K ( t 



-7-1 



(20) 



with < 7 < 1. Here k is a dimensionless constant and 
To is a characteristic mesoscopic time (but not a first 
moment). The small-it behavior of the Laplace transform 
of the waiting time distribution reads 



ip(u) ~ 1 - Au 1 , 



(21) 



with A = 7~ 1 f«r(l — J j i D . 

We take the jump distance distribution to be "nor- 
mal," that is, it has finite moments. In this case, we are 
interested in small values of q = |q|, for which one can 
expand the Fourier transform of the jump displacement 
distribution and retain only the first two terms, 



x(q) ~ i 



a 2 q 2 



where the second moment 



a 2 = J drr 2 x(r) 



(22) 



(23) 



is assumed to be finite. 

Substitution of the expansions ([2Tj) and (|22"|) into 
Eq. (fTT)|) and neglect of a term of 0(itq 2 ) (which is unim- 
portant in the asymptotic regime of small wave vectors 
and low frequencies) leaves us with 



urj(q., u) 



1—1° " ~l \ 



¥>(ro,0) 

Laplace and Fourier inversion then yield 
0»?(r,f|r O) O) 



dt 



(24) 



(25) 



where we have introduced the anomalous diffusion coef- 
ficient 



K 7 = ^. (26) 

The integrodifferential operator T> t 7 acting on y(t) is 
defined as the inverse Laplace transform of u 1_7 y(w), 



(27) 



and is closely related to the Riemann-Liouville operator 



qD} 7 defined by 



1 d 



dt' 



(28) 



In fact, both operators V\ 7 and qD] 7 are the same 
when applied to sufficiently regular functions f(t) as de- 
termined by the condition lim t ^o Jq d,T(t — t) 7 _1 /(t) = 
(see Pgs. 384 and 118 in [l^). This condition is satisfied 
for all situations of interest here. 

In terms of the probability density w, Eq. (|2"5j) explic- 
itly yields the fractional reaction-subdiffusion equation 



dw{r,t\r o ,0) 
dt 



= <p(r,t)K 1 oVl^Vl 



1 



I — 



u;(r,t|r o ,0) 



V>(r,t) 



w(r,t\r o ,0). 



(29) 



Equation (|2^|) is the first mesoscopic highlight of our pa- 
per in that all further results for subdiffusion are ob- 
tained as special cases. It is perhaps the most general 
fractional subdiffusion equation associated with a CTRW 
with evanescence for a single-species system obtained to 
date. It is not as general as one might hope because of 
the rather stringent condition (fT0|) . but to the best of 
our knowledge it does include the equations currently in 
the literature. The results of Sokolov et al. Q and the 
explicit results in Henry et al. Q are recovered if we set 
ip(r, t) cx exp(— kt). 



Somewhat more elaborate is the connection with the 
recent work of Fedotov [3, which in turn recovers the 
more general results in 8]. In the language of Fedotov 
translated to our work, his choice corresponds to the par- 
ticular selection (in one dimension) 



ip(x, t) = exp < / r[p(x,t')]dt' 



(30) 



Here in Fedotov's language the chemical reaction respon- 
sible for the evanescence is assumed to follow the law of 
mass action so that the reaction term is of the form r(p)p, 
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that is, 



r[p{x,t)}. 



(31) 



Reaction 



In turn, p(x, t) is the density of (surviving) particles at 
point x at time t, which is related to our probability 
density via an integration over all initial positions, 



p(x,t) = / dxotv(x, t; xq, 0). 



(32) 



The time r in Eq. (|30p is an arbitrary reference time 
that can be chosen to be zero. Specifically, the choice 
r = 0, substitution of Eq. (|30| into our general equa- 
tion (|29|) , and subsequent integration over the initial con- 
dition leads exactly to Eq. (19) in Ref. 0, which is thus 
again a special case of our more general formalism. The 
result is especially noteworthy because the probability 
of spontaneous death is not set a priori but depends on 
the changing density itself, This leads to an interesting 
complex nonlinear problem. 

As a final note in this section we mention that under 
the same conditions that led to our subdiffusive fractional 
equation (|29p with a Riemann-Liouville operator, we can 
arrive at an equivalent fractional subdiffusion equation 
of Caputo form. It turns out to be 



di w(r,t|ro,0) 
<9iT tp(r, t) 



^M|ro,0) 



(33) 



This is in general not our preferred choice because of the 
difficulties in carrying out the Caputo fractional deriva- 
tive of a product. 



B. Fractional superdiffusive equation with 
evanescence 

In the previous subsection we dealt with a waiting time 
distribution with long tails together with a "normal" dis- 
tribution of displacements, and the end result was a frac- 
tional subdiffusion equation. In this subsection we con- 
sider a "normal" waiting time distribution along with a 
distribution of displacements that has long tails. This 
will lead to a superdiffusive equation. 

We thus consider a waiting time distribution with a 
finite mean waiting time r between steps, so that its 



Laplace transform at small argument (corresponding to 
long times) behaves as 



1 — TU. 



(34) 



For the jump length distribution we assume an inverse 
power law behavior, 



x( r ) 



(35) 



with 1 < fi < 2 and r = |r|. Its variance diverges, and 
its Fourier transform is 

*(q) ~ 1 - a*?", (36) 
where q = |q|. The steps to follow are now straightfor- 
ward. Again, it turns out to be convenient to work with 
Eq. (|17[) . Substitution of the above expansions into this 
equation, retention of leading terms, and some simple 
algebra leads to 



uf/{q,u) - 



¥>(ro,0) 

or, inverting the time Laplace transfrom, 
d 

— f)(q,t) = -—^(q,*). 



— ^(q,«), (37) 

T 



(38) 



Let d^/dr^ be the operator defined by the following 
Fourier transform property, 



We can thus write 

at 



(39) 



= K— 77(r,i|r ,0), (40) 



where we have introduced the anomalous diffusion coef- 
ficient 



T 



(41) 



For a one-dimensional system d^/dr^ = d^/dx^" is the 
Riesz operator [To|, HH, Ell • 

Undoing the variable transformation (|16p . we finally 
arrive at the fractional reaction-superdiffusion equation, 



dw(r,t|r ,0) 
dt 



= <P(T,t)K 



1 



w(r,i|r o ,0) 



<P(T,t) 



w(r,t\r ,0). 



(42) 



This is the second important mesoscopic result of our fractional equation starting from a CTRW. 
paper, namely, the derivation of a reaction-superdiffusion 
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C. Bifractional equation with evanescence 

Finally, in this subsection we combine subdiffusion and 
superdiffusion in that we choose a waiting time distribu- 
tion that lacks finite moments (and thus leads to subd- 
iffusion by itself) with a jump distance distribution that 
also lacks moments (and thus leads to superdiffusion by 
itself). Our methodology directly lends itself to this com- 
bination. 

We choose the waiting time distribution of Eq. ([20]) 
whose Laplace transform is given in Eq. (f2Tj) . and the 



J 



jump distribution Eq. f|35[) whose Fourier transform is 
given in Eq. (|36|) . The steps to follow are now essentially 
the same as in the previous sections, with appropriate 
care given to the retention of the leading contributions. 
After some algebra we find 

«£^ = * oP ;-|>, t K,0). (43) 



Undoing the variable transformation (|16[) . we finally ob- 
tain the fractional reaction-sub/super-diffusion equation, 



dw(r, t\r , 0) 

di 



= <p(r, t) K V\- 



d" 



1 



dr^ ip(r, t) 
I — 



w(r,i|r o ,0) 



y(r,t) 



w(r,t\r ,0). 



(44) 



This is our third mesoscopic result and is unique in that it 
combines both subdiffusion and superdiffusion in a single 
equation. 



V. SUMMARY AND OUTLOOK 

In this paper we have approached the problem of de- 
scribing the evolution equation of particles that move in 
a medium in which they can also die in a medium as they 
move. The model is based on a CTRW description of the 
motion of the particles. The motion may be anomalous 
(subdiffusive or superdiffusive) and the particles may die 
at a rate that can depend on position as well as time. We 
are able to capture the models that have been explicitly 
considered in the literature, e.g., the space- independent 
exponential evanescence model of Refs. [3j and [8[, but 
our model can also capture complicated position depen- 
dences of the evanescent behavior such as that of the 



model of Fedotov L| that render the problem nonlinear. 
We confirm in a more general way than had been es- 
tablished previously the known result that the interplay 
of the (anomalous) motion and the evanescence is quite 
complex and that in general it can not be represented 
as the sum of two processes the way it can in normal 
reaction-diffusion scenarios. Having said this, we note 
as an aside that in certain cases (e.g., that of evanes- 
cence at a constant rate r = k) it is possible to reduce 
the reaction-subdiffusion problem to a pure subdiffusion 
problem by a proper transformation [14j . much in the 
spirit of Danckwerts' solution for the problem of classical 
diffusion with a linear reaction [15J. 

We started by constructing an integral equation for the 
probability density of finding a (surviving) particle at a 
location r at time t given that it stepped on location ro at 
time t = 0. To proceed from this CTRW-based equation 
to the fractional equations, we found it necessary to make 
some specific assumptions about the form of the rate of 



evanescence. The kernel of the integral equation under 
these conditions is simply related to the single-step jump 
probability density of the underlying CTRW. In spite of 
the constraints, our models include as special cases all the 
explicit models that have been presented in the literature. 
All those, case by case have been presented with 

very specific physical contexts in mind. 

Once we have arrived at an integral equation, the 
derivation of various fractional diffusion equations relies 
on fairly standard procedures dictated by the form of the 
single-step probability properties, except that we have 
added evanescence to the picture and are thus able to see 
the complex interplay of motion and evanescence, at least 
under our assumptions. Eventually we hope to be able to 
relax some of our more stringent assumptions. We also 
hope to be able to include other particle loss mechanisms 
such as bimolecular reactions in our scheme, possibly at 
the expense of introducing some kind of mean-field as- 
sumption to deal with the complexity arising from effects 
of cooperativity. The existing formalism can be adapted 
to some situations where instead of evanescence we have 
particle sources, or perhaps sources and sinks simulta- 
neously. However, in some cases particle sources require 
special scrutiny. For example, special care is needed when 
dealing with particles that give rise to offspring because 
one must specify the rules surrounding the location and 
time of creation of new particles, especially when deal- 
ing with jump and waiting time distributions that have 
long tails. These are all plans for future work. Our most 
immediate plans are to apply our results to the problem 
of the survival probability of a target that is surrounded 
by a d-dimensional sea of evanescent traps whose motion 
may be subdiffusive if the waiting times for motion are 
too long, or superdiffusive if the jumps are sufficiently 
long, or a mixture of both. 

Finally, we end with an interesting observation that 
does not appear obvious. One might, instead of Eq. (jT4j) , 
have been motivated to write the integral equation 
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u>(r,i|r ,0) = J dr' J dt' *(r' - r 0) i')H(r', r ;t', 0)u;(r, t\r', t') + $(i)H(r 0) r ; t, 0)6(r - r ). 



(45) 



r 



Again in the language of Hughes [l3(, this corresponds 
to a "partition over the first step." What is interesting is 
that we are not able to arrive at any reasonable fractional 
diffusion equation starting from this integral equation, in 
any case not by the methods followed in this paper for 
Eq. ITU) even though one might have expected a certain 
symmetry to the situation. This, too, is a question to be 
explored further. 
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